Observational constraints on late-time A(£) cosmology 
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The cosmological constant A, i.e., the energy density stored in the true vacuum state of all 
existing fields in the Universe, is the simplest and the most natural possibility to describe the current 
cosmic acceleration. However, despite its observational successes, such a possibility exacerbates the 
well known A problem, requiring a natural explanation for its small, but nonzero, value. In this 
paper we study cosmological consequences of a scenario driven by a varying cosmological term, in 
which the vacuum energy density decays linearly with the Hubble parameter, A oc H. We test 
the viability of this scenario and study a possible way to distinguish it from the current standard 
cosmological model by using recent observations of type la supernova (Supernova Legacy Survey 
Collaboration) , measurements of the baryonic acoustic oscillation from the Sloan Digital Sky Survey 
and the position of the first peak of the cosmic microwave background angular spectrum from the 
three-year Wilkinson Microwave Anisotropy Probe. 



I. INTRODUCTION 

The nature of the mechanism behind the current cos- 
mic acceleration constitutes a major problem nowadays 
in Cosmology [1]. Even though almost all observational 
data available so far are in good agreement with the sim- 
plest possibility, i.e., a vacuum energy plus cold dark 
matter (ACDM) scenario, it is becoming rather consen- 
sual that in order to better understand the nature of the 
dark components of matter and energy, one must also 
consider more complex models as, for instance, scenarios 
with interaction in the dark sector Q. 

In this regard, the simplest examples of interacting 
dark matter/dark energy models are scenarios with vac- 
uum decay [A(t)CDM]. In reality, A(t)CDM cosmologies 
constitute the special case in which the ratio of the dark 
energy pressure to its energy density, uj, is exactly —1 
@j W- This kind of model may be based on the idea that 
dark energy is the manifestation of vacuum quantum fluc- 
tuations in the curved space-time, after a renormalization 
in which the divergent vacuum contribution in the flat 
space-time is subtracted. The resulting effective vacuum 
energy density will depend on the space-time curvature, 
decaying from high initial values to smaller ones as the 
universe expands [5J. As a result of conservation of to- 
tal energy, implied by Bianchi identities, the variation of 
vacuum density leads either to particle production or to 
an increasing in the mass of dark matter particles, two 
general features of decaying vacuum or, more generally, 
of interacting dark energy models @ . 

Naturally, the precise law of vacuum density variation 
depends on a suitable derivation of the vacuum contri- 
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bution in the curved background, which is in general a 
difficult task. In this regard, however, a viable possi- 
bility is initially to consider a de Sitter space-time and 
estimate the renormalized vacuum contribution with help 
of thermodynamic reasonings as, e.g., those in line with 
the holographic conjecture [7[. The resulting ansatz is 
A w (H + to) 4 — to 4 , where H is the Hubble parameter 
and to is a cutoff imposed to regularize the vacuum con- 
tribution in flat space-time (the next step is to consider 
a quasi-de Sitter background, with a slowly decreasing 
H. In this case, the above ansatz may be considered a 
good approximation, but with the vacuum density decay- 
ing with time.). In the early-time limit, with H » to, 
we have A « H 4 . By using this scaling law for the vac- 
uum density and introducing a relativistic matter com- 
ponent, some of us have obtained from the Einstein equa- 
tions an interesting solution with the following features 
[gj. Firstly, the Universe undergoes an empty, quasi-de 
Sitter phase, with H « 1, which in the asymptotic limit 
of infinite past tends to de Sitter solution with H = 1 
(in Planck units). However, at a given time, the vacuum 
undertakes a fast phase transition, with A decreasing to 
nearly zero in a few Planck times, producing a consid- 
erable amount of radiation. In this sense, this can be 
understood as a non-singular inflationary scenario, with 
a semi-eternal quasi-de Sitter phase originating a radia- 
tion dominated universe. 

The subsequent evolution of vacuum energy essentially 
depends on the masses of the produced particles. If we 
apply the energy-time uncertainty relation to the pro- 
cess of matter production, we will conclude that mas- 
sive particles can be produced only at very early times, 
when H is very high. Therefore, baryonic particles and 
massive dark matter (as supersymmetric particles and 
axions) stopped being produced before the time of elec- 
troweak phase transition. On the other hand, the late- 
time production of photons and massless neutrinos must 
be forbidden by some selection rule, otherwise the Uni- 
verse would be completely different from the one observed 
today. Thus, if no other particle is produced, the vac- 
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uum density stops decaying at very early times, and fol- 
iate times we have a standard universe, with the pres- 
ence of a genuine cosmological constant (some thermody- 
namic considerations, in line with the holographic prin- 
ciple, permit to infer the value of this constant, leading 
to A « to 6 [li]). 

In order to have a decaying vacuum density at the 
present time, the produced dark particles should have 
a mass as small as 10 _65 g 1 . Here, we have considered 
the possibility of a late-time decaying A, and compared 
the consequent cosmological scenario with the constraints 
imposed by current observations [Tol.[T]|. In such a limit, 
H « m and A w m 3 H, leading again to A w m 6 in the 
de Sitter limit. From a qualitative point of view, we have 
found no important difference between this A(i) scenario 
and the flat ACDM model 10]. After the phase tran- 
sition described above the Universe enters a radiation- 
dominated phase, followed by a matter-dominated era 
long enough for structure formation, which tends asymp- 
totically to a de Sitter universe with vacuum dominating 
again. The only important novelty, related to matter 
production, is a late-time suppression of the density con- 
trast of matter, which may constitute a potential solution 
to the cosmic coincidence problem [12j. Moreover, the 
analysis of the redshift-distance relation for supernovae 
of type la, particularly with the Supernova Legacy Sur- 
vey (SNLS) data set, has shown a good fit, with present 
values of H and the relative density of matter in accor- 
dance with other observations [111 ]. 

In this paper, we go a little further in our investiga- 
tion and study new observational consequences of the 
A(t)CDM scenario described above. We use to this end 
distance measurements from type la supernovae (SNe 
la), measurements of the baryonic acoustic oscillations 
(BAO) and the position of the first peak of the cosmic 
microwave background (CMB). We show that, besides 
the interesting cosmic history of this class of A(t)CDM 
models, a conventional, spatially flat ACDM model is 
only slightly favored over them by the current observa- 
tional data. 



II. THE MODEL 

In a spatially flat FLRW space-time the Friedmann and 
the conservation equations can be written, respectively, 
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as 



and 



p T = 3H 2 



pr + 3H(p T + pt) = 



(1) 



(2) 



1 Some authors associate this mass to the graviton in a de Sitter 
background with H and A of the orders currently observed [tj . 

2 We work in units where Mpi anc i< = (87TC?) -1 / 2 = c = 1. 



where px and pt are the total energy density and pres- 
sure, respectively. If we consider that the cosmic fluid is 
composed of matter with energy density p m and pressure 
p m , and of a time-dependent vacuum term with energy 
density p\ = A and pressure pa = —A, we obtain 



3H( Pr , 



Pn 



-A 



(3) 



which shows that matter is not independently conserved, 
with the decaying vacuum playing the role of a matter 
source. Throughout our analysis we assume that baryons 
are independently conserved at late times, being not pro- 
duced at the expenses of the decaying vacuum. This 
amounts to saying that we will postulate, in addition 
to Eq. ([21), a conservation equation for baryons, i.e., 
pb + 3H(pb + pi,) = 0, where pb and pb refer to baryon 
density and pressure. 

From the above equations and considering our late- 
time ansatz A = aH (a is a positive constant of the 
order of to 3 ), the evolution equation reads 



2H + 3H 2 -<rH = 



(4) 



Now, with the conditions H > and p m > 0, the 
integration of the above equation leads to the following 
expression for the scale factor, 



a{t) = C [exp (at/2) - lp 



(5) 



where C is the first integration constant and the second 
one was chosen so that a = at t = 0. From these equa- 
tions, it is straightforward to show that at early times 
the above expression reduces to the Einstein-de Sitter 
solution whereas at late times it tends to the de Sitter 
universe. 

As stated earlier, A = aH and, therefore, p m — 3H 2 — 
aH. By using solution ([5]), one can also show that 



7 2 C 3 a 2 C^ 2 
3a 3 + 3a 3 / 2 



and 



A = — 
3 



a 2 . a 2 C 3 ' 2 
3a 3 / 2 ' 



(G) 



(7) 



The above expressions can be easily understood as fol- 
lows. The first terms in the r.h.s. are the expected scaling 
laws for matter density and the cosmological constant in 
the case of a non-decaying vacuum while the second ones 
are related to the time variation of the vacuum density 
and the concomitant matter production. As expected, 
at early times matter dominates, with its density scaling 



with 



and the matter production process is negligi- 



ble. On the other hand, at late times the vacuum term 
dominates, as should be in a de Sitter universe. 

From Eqs. |T]), ^ and (O, the evolution of the Hubble 
parameter as a function of the redshift can be written as 



H{z) = H l-Q m + Q m (l + z) 



.3/2 



(8) 
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where Hq and Q m are, respectively, the present values of 
the Hubble parameter and of the relative energy density 
of matter. Note that, due to the matter production, this 
expression is rather different from that found in the con- 
text of the standard ACDM case. In particular, if A = 
and fl m = 1; we obtain H{z) = Hq(1 + z) 3 / 2 , leading to 
p m = 3H 2 (1 + z) 3 , as expected for the Einstein-de Sitter 
scenario. 



TV = 115. As discussed in Ref. [TO], the best-fit 
values for this analysis is obtained for h — 0.70 and Q m — 
0.33, with reduced Xr = 1-01. At 95% of confidence 
level, we also find 0.69 < h < 0.71 and 0.28 < Q m < 
0.37. The 2-dimensional contours shown in the Figures 
are obtained from the traditional frequentist confidence 
intervals (based on the A% 2 approach and assuming that 
errors are normally distributed). 



III. OBSERVATIONAL CONSTRAINTS 



B. Baryonic acoustic oscillations 



Extending and updating previous results [ll( , we study 
in this Section some observational consequences of the 
class of A(t)CDM scenarios discussed above. Note that, 
similarly to the standard ACDM case, this class of models 
has only two independent parameters, H and fi m [see 
Eq. ©]. The best-fit values for these quantities will be 
determined on the basis of a statistical analysis of recent 
type la supernovae (SNe la) measurements, as given by 
the SNLS Collaboration [l3j , the distance to the baryonic 
acoustic oscillations (BAO) from the Sloan Digital Sky 
Survey (SDSS) [lij], and the position of the first peak 
in the spectrum of anisotropies of CMB radiation from 
the three-year Wilkinson Microwave Anisotropy Probe 
(WMAP) [15| (for more details on the statistical analysis 
discussed below we refer the reader to Ref. [l6| ) . 



A. SNe la observations 

The predicted distance modulus for a supernova at red- 
shift z, given a set of parameters s, is 



fj, p (z\s) = m — M = b\ogd L + 25, 



(9) 



where m and M are, respectively, the apparent and ab- 
solute magnitudes, the complete set of parameters is 
s = (Hq, fi m ), and stands for the luminosity distance 
(in units of megaparsecs), 



d L = c(l + z) 



dx 



x , x 2 H(x;s)' 



(10) 



with x' = ?jp- = (1 + z) 1 being a convenient integration 
variable, and H(x; s) the expression given by Eq. (JS}. 

We estimated the best fit to the set of parameters s by 
using a x 2 statistics, with 



A' 



[/4(*|s)-/4(*)] : 



(11) 



where [j? (z\s) is given by Eq. ([9]), p l {z) is the extinction 
corrected distance modulus for a given SNe la at Zj, and 
°i = a (^B) + °int> where °f MB ) is the variance in the 
individual observations and a 2 nt stands for the intrinsic 
dispersion for each SNe la absolute magnitude. Since 
we use in our analysis the SNLS collaboration sample 



The use of BAO to test dark energy models is usually 
made by means of the parameter A, i.e., [l4j 



A = A 



(12) 



where z — 0.35 is the typical redshift of the sample and 
Dy is the dilation scale, defined as 



D v (z) = 



D M {zf 



H(z) 



1/3 



with the comoving distance Dm given by 



D M (z) = c 



dz' 
o W) 



(13) 



(14) 



An important aspect worth emphasizing at this point 
is that the value of A is obtained from the data in the 
context of a ACDM model, and can be considered a good 
approximation only for some class of dark energy models 
|17| . In particular, two conditions are implicitly assumed 
to be valid [T3|. First, the evolution of matter density 
perturbations during the matter-dominated era must be 
similar to the ACDM case, at least until the character- 
istic redshift z = 0.35. Second, the comoving distance 
to the horizon at the time of equilibrium between mat- 
ter and radiation must scale with (f2 m i?g) • However, 
none of the above conditions are satisfied in the present 
model because of the matter production associated to the 
vacuum decay. As will be shown in a forthcoming pub- 
lication [ijj , if matter is homogeneously produced there 
is a suppression in its density contrast for z < 5, that is, 
after the period of galaxy formation (this will eventually 
imply a higher value of £l m in order to fit the observed 
mass power spectrum). 

On the other hand, as radiation is independently con- 
served, its relative energy density for z >> 1 is given 
by f2 r z 4 , where tt r is its present value. With the help 
of |8]) we can see that, in the same limit z >> 1, 
the relative density of matter is Fl^z 3 (with the ex- 
tra factor f2 m being due to the matter production be- 
tween t(z) and the present time). By equating the 
two densities, we obtain the redshift of equilibrium be- 
tween matter and radiation, given by z rm — fl 2 n /il r . 
Therefore, after including conserved radiation into ([5]) 
(see equation (|17l) below) and taking z >> 1, we have 
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fH(zrm) ~ Cy/ p r /Q(fl m Ho) 2 , where p r is the present 
value of the radiation density (while in the ACDM case 

we would obtain ru{z rm ) « CyJ p r / '6(f2 m ffo) _1 > as stated 
above). Thus, one can see that the parameter A is not 
appropriate to test the model, and we will use instead the 
dilation scale Dy, which is weakly sensitive to the cos- 
mological evolution before z — 0.35. By combining our 
function H(z), given by |5]), into (fT5|) - p^|) we can find 
the region in the fi m — Ho plane which gives the observed 
value D v = (1370 ±64) Mpc (la) [H|. The BAO bands 
in the fl m — Hq parametric space are shown in Fig. 1. 



parameters at high redshifts, due to the process of matter 
production. Therefore, in order to perform a correct test, 
we have to explicitly calculate the acoustic scale in the 
model and then compare with the measured position of 
the peak. 

The acoustic scale, defined as the ratio between the 
comoving distance to the surface of last scattering and 
the radius of the acoustic horizon at that time, is then 
given by 



C. The first peak of CMB 

The two tests we have previously described depend on 
the physics at low and intermediate redshifts (until z ~ 1) 
and will lead, as we will see, to very similar results. As a 
complementary test, involving high-z measurements, let 
us consider the position of the first peak in the spectrum 
of CMB anisotropies. Since in the present model there 
is no production of baryonic matter or radiation and the 
spatial curvature is null, we expect that a correct posi- 
tion of the first peak is enough to guarantee a spectral 
profile similar to the ACDM one, provided the spectrum 
of primordial fluctuations is the same. 

In the context of a large class of dark energy models 
this test is performed by comparing the predicted shift 
parameter with the ACDM value. However, this is only 
valid if the acoustic horizon at the time of last scattering 
is the same [l8|. This is not true in the present model 
because, as we have already shown, for the same values of 
Hq and Q m we have different expressions for cosmological 



r z is dz 
, _ 71 JO H(z) 
A ~ r°° c s dz 
Jz u c H(z) 



where zi s is the redshift of last scattering, and 



/„ , 9 n b \ 

c s = c 3 + - — — 
\ 4 fl 7 z J 



is the sound velocity. Here, fib and 7 are, respectively, 
the present relative energy densities of baryons and pho- 
tons. 

The function H(z) to be used in Eq. (|15[) must now 
include radiation. As this component is independently 
conserved, scaling with a -4 , the appropriate generaliza- 
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FIG. 2: Left) Confidence regions (68.3%, 95.4% and 99.7%) in the Q m — h plane for a joint analysis involving SNe la, BAO 
and CMB data. As indicated, the dashed lines stand for the SNe la results of Ref. [llj]. In both panels, the largest contours 
represent the joint analysis of SNe la plus BAO measurements whereas the smallest ones arises when the CMB data are included 
in the analysis. Right) The same for a spatially flat ACDM model. 



tion of Eq. {SI is given by 3 



H(z) 
Ho 



i -n m + n m (i + z) 



3/2 



n r (i 



1/2 



(17) 

Therefore, apart from our two free parameters Hq and 
Q m , in order to determine the acoustic scale we need the 
present values of the energy densities of radiation, pho- 
tons and baryons, as well as an expression for z\ a . Since 
radiation and baryons are independently conserved, and 
we want to preserve the spectrum profile as well as the nu- 
cleosynthesis constraints, we will take for these densities 
the same values obtained from CMB observations in the 
context of the ACDM model H 7 /i 2 « 2.45 x 1(T 5 , 
n r h 2 w 4.1 x 1CT 5 and Q b h 2 « 0.02. 

Concerning zi s , its value is not in principle the same 
as in the ACDM case. To determine its value for a given 
pair (n m ,ifo) we will proceed as follows. First, we obtain 
the redshift of last scattering z ; * s in the ACDM model by 



3 The inclusion of a conserved component of radiation changes the 
dynamics and, consequently, the evolution of A(jz) and p m (z). 
Thus, rigorously speaking the generalization of H(z) would re- 
quire a reanalysis of the dynamics. Nevertheless, as £7 r fa 10~ 4 , 
when the decaying vacuum and matter production begin to be 
important, the radiation term is negligible. For this reason, Eq. 
1171 can be considered a good approximation. Indeed, a numeri- 
cal analysis in the range < z < 10 4 showed a difference between 
Eq. 117> and the exact H(z) as small as 0.01%. 



means of the current fitting formula [l8j, and then we 
impose that the optical depth has to be the same in both 
models, that is, 



T(z) dz 
H{z) l + z 



r(z) dz 
H*(z)l + z' 



(18) 



where T is the rate of photon scattering and H* is the 
Hubble parameter in the ACDM context. For the current 
intervals of Hq and Q m , the relative differences between 
zi s and z* s are typically as small as 1%. 

For a scale invariant ACDM model with spectral index 
n = 1, the position of the first peak after including the 
effect of plasma driving is given by the fitting expression 

M 



h = l A (l-6 1 ), 



where 



Si = 0.267 



f— V 
V0.3/ 



(19) 



(20) 



with r = p-y(zi s ) I 'pm(zis) ■ Since the plasma driving de- 
pends essentially on pre-recombination physics, we will 
therefore assume that the above fitting formulae are a 
good approximation to the present A(t) model. In our 
case, the parameter r is given by 



fil 



-Zis 



(21) 
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Finally, from the above results we can determine the 
region of the f2 m -iJ plane for which the first peak has the 
position currently observed by WMAP, i.e., l\ — 220.8 ± 
0.7 (la) As in the case of BAO, the CMB bands in 
the O m — Hq parametric space is shown in Fig. 1. 



D. Results and discussions 

The superposition of the confidence regions and bands 
corresponding to our analysis of SNe la, BAO and CMB 
observations is shown in Figure 1. For the sake of com- 
parison, the right panel shows the same analysis for the 
standard ACDM model. From these analysis, we note 
that, although not very restrictive and parallel to the 
SNe la contours, the BAO bands can be statistically 
combined with SNe la data providing constraints on the 
Vt m - H plane. At 2a level, we find h = 0.70 ± 0.01 and 
Q m = 0.32 ±0.05, with reduced xf - 1-00. The inclusion 
of CMB data into our analysis in turn makes the com- 
plete joint analysis more restrictive but also shows that 
the concordance is as good as in the standard ACDM 
case, and that the A(t)CDM model discussed here can- 
not be ruled out. It can be anticipated from Figure 1 
that the data will prefer higher values of the matter den- 
sity parameter, while the Hubble parameter is expected 
to be slightly smaller than the current accepted values. 
In order to confirm this qualitative discussion, Figure 2 
shows the results of our joint statistical analysis (SNe la 
+ BAO + CMB). At 2a level, we find h = 0.69±0.01 and 
Q m = 0.36 ± 0.01, with reduced xl = 1-01. For the sake 
of comparison, the same analysis for the ACDM case is 
also shown in the right panel of Fig. 2. 

Finally, with the above best-fit values for the param- 
eters, we can calculate the total ex pan ding age for this 
class of A(t)CDM models, given by [10L 1111] 

io = ^ 1 ^r-15.0Gyr, (22) 

as well as the redshift of transition from a decelerated 
expansion to the current accelerating phase, i.e., zt ~ 
1.3. Note that both values are slightly higher than, but 
of the same order of, the standard ones. 



IV. FINAL REMARKS 

By using the most recent cosmological observations, 
we have discussed the observational viability of a class 
of A(t)CDM scenarios in which A oc H, as well as a 
possible way to distinguish it from the standard ACDM 
model in what concerns the general characteristics of the 
predicted cosmic evolution. As discussed earlier, these 
A(t)CDM models have some interesting features as, e.g., 
the association of dark energy with vacuum fluctuations, 
the circumvention of the cosmological constant problem 
by subtracting the flat space-time contribution from the 



curved space-time vacuum density, and the possible (but 
not necessary) link between dark matter and massive 
gravitons. 



We have presented some quantitative results which 
clearly show that, even in the current stage of the Uni- 
verse evolution, our decaying vacuum scenario is very 
similar to the standard one. We have statistically tested 
the viability of the model by using recent SNe la observa- 
tions and measurements of the BAO and the first peak of 
the CMB spectrum. At 95.4% c.L, a joint analysis involv- 
ing SNe la + BAO provides the intervals 0.69 < h < 0.71 
and 0.28 < ft m < 0.37, which are in good agreement 
with the values of the Hubble and the matter de nsity 
parameters obtained from independent analysis |2ll. l22j|. 
When the position of the first peak of CMB anisotropics 
is included in the analysis, the best-fit value for the rel- 
ative matter density is increased, fl m ~ 0.36. This re- 
sult cannot rule out the model, and it may be indicating 
that, besides the interesting cosmic history of this class 
of A(t)CDM models, a conventional, spatially flat ACDM 
model is only slightly favored over them by the current 
observational data. Still on the best-fit for fi m , we note 
that such a higher matter density is something to be more 
investigated by means of other cosmological or dynam- 
ical tests as, e.g., the predicted mass power spectrum 
in the context of the model. As we have discussed ear- 
lier, a higher matter density is necessary to compensate 
the late-time suppression of the density contrast owing 
to matter production (T^ |. Another possibility will be 
provided by future supernovae observations, since the 
present model starts to diverge from the standard one 
for higher redshifts [TTI ] . 



Finally, we should also emphasize two aspects to be 
considered before a definite conclusion about the com- 
parison between the present model and the flat standard 
scenario. First, that in our study of CMB we have used 
parameter values and fitting formulae that are strictly 
correct for the ACDM case, particularly the expression 
giving the position of the first peak for a given acoustic 
scale [equation (flU)) ]. In spite of our reasons to consider 
that use as a good approximation, it can lead to bias 
in our results. If that is the case, only a more com- 
plete analysis of CMB would rule out or corroborate the 
model. The second aspect is of theoretical character. As 
discussed in the Introduction, our ansatz for the variation 
of A, if genuine, is a good approximation only for quasi- 
de Sitter backgrounds. Naturally, our Universe, although 
dominated by a cosmological term, is far away from the 
asymptotic de Sitter state, with matter still giving an 
important contribution to the cosmic fluid. Only a more 
profound theoretical study, on the basis of quantum field 
theories in the expanding background, could establish the 
degree of applicability and limits of that approximation. 
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